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An algorithm for optimizing the MERA tensor network in an infinite system is presented. Using 
this technique we compute the critical exponents of Ising and XXZ model. 
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I. INTRODUCTION 

Critical phenomena are ubiquitous in science ranging 
from condensed matter to biological or economic systems. 
They are associated to scale invariance, a diverging cor- 
relation lenght, and the various correlation functions de- 
cay as power law. A key issue in the study of of criti- 
cal systems is the computation of the critical exponents, 
i.e. the exponents which govern the decay of the corre- 
lations. The power of methods like the renormalization 
group stems from its capability to compute (analitically 
or numerically) the critical exponents. 

In the field of quantum many-body systems important 
progresses have been made recently in one-dimension 
where a number of tensor network approximations [J, [U, 
& 0, IE B S Hi °f the many-body wave- function allowing 
an efficient computation of all the relevant observables. 
The Multi-scale Entanglement Renormalization Ansatz 
(MERA) introduced by Vidal [8( is particularly appealing 
since embeds scale invariance typical of critical systems. 
A method to relate the critical exponents in the MERA 
tensor network has been very recently put forward by 
some of us [9j by relating them to the eigenvalues of the 
MERA transfer matrix. 

In the present work we want to apply this method for 
the computation of the critical exponents of two corner- 
stone models in quantum statistical mechanics, the Ising 
and the XXZ model. To this end we provide an algorithm 
for optimizing the MERA tensor network in an infinite 
system and then, by employing the results given in [9( we 
compute the critical exponents. 

The paper is organized as follows. In the next section 
we briefly review the optimization of the MERA tensor 
network and how to compute the critical exponent. We 
then proceed with the comparison with the Ising and 
XXZ model defined in Section Hill In Section ITVl we dis- 
cuss the results and compare them with the exact known 
values. The last Section is devoted to the conclusions. 



II. MERA OPTIMIZATION 

Consider a 1-D translational invariant many-body 
quantum system at criticality composed of L = 2 i+1 



sites, each of them described by a local Hamiltonian and 
some nearest neighbor interactions. The global Hamil- 
tonian can then be expressed as H = ^2 k h k , where 
k = 1, • • • , L with L+ 1 = 1 for periodic boundary condi- 
tions and where h k is acting on the three consecutive sites 
k — 1, k, and k + 1 (the generalization to longer range in- 
teraction is straightforward) . In the thermodynamic limit 
of infinitely many sites the ground state energy per site 
of the system can thus be computed as 
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where |\&) are joint states of the many-body system, and 
where the last minimization is performed over the sets 
By := {p\, ■ ■ ■ , pl} whose elements can be obtained as 
reduced density matrices of some global pure state |\&) of 
the system (specifically p k is the reduced density matrix 
of that is associated with the sites k — 1, k and k + 
1). Owing to translational invariance all the addends 
of the last term of Eq. ((!]) are equal. Thus without the 
constraint p k € By the minimization would be trivial and 
Eq will coincide with the the ground energy level of the 3- 
sites Hamiltonian hk- Taking into account the condition 
Pk E By is what makes the calculation of Eq a hard 
problem to solve: it requires us to minimize the energy 
of hk by properly embedding p k in the "environment" 
formed by the remaining sites of the many-body system. 

A solution can be found by adopting the so 
called Multi-scale Entanglement Renormalization Ansatz 
(MERA) [8| , which, similarly to the matrix product state 
decomposition [l|, 0, H, 0, [g], assumes a specific repre- 
sentation of the many-body state |\&). In particular ac- 
cording to this ansatz the wave function of an L sites 
many-body system is expressed in terms of 0(Llog 2 L) 
tensors of two different species (i.e. the type-( \ ) disen- 
tangler tensors x an d type-( \ \ isometry tensors A) which 
are connected to form a multi-layer structure which ad- 
mits efficient contraction rules (see Fig. [1] — we refer 
the reader to Ref. [8( for details). Specifically, following 
Ref. Q we will restrict the minimization (Q]) to states |\&) 
which can be expressed (or at least approximated) by ho- 
mogenous MERAs, in which all the disentangler \ an d 
all the isometry A entering the tensor network are identi- 
cal. Under these conditions it has been shown Q that the 
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thermodynamic limit L — > oo of |\&) can be characterized 
in terms of a transfer super-operator (the QuMERA map 
<&) that is determined by the x s an d the As that forms 
the MERA [To| ■ In particular the reduced density matrix 
Pk of | 'J) can be computed as the (unique) eigenvector p* k 
of the MERA transfer super-operator (QuMERA map) 
$ associated with the unitary eigenvalue, i.e. 

m)=pt- (2) 

This means that Eq. |T]) can now be expressed as 

E G = min$ Tz\p* k h k ] , (3) 

where the minimization is performed over the set of all 
possible QuMERA channels $ - the equivalence being 
guaranteed by the stationary condition ([2]) . Equation ([3]) 
achieves two fundamental goals: i) it allows us to di- 
rectly address the thermodynamic limit, ii) it guaran- 
tee the possibility of reconstructing a many-body joint 
state | \Jj) associated with the three sites local density 
operator p* k (this is the MERA state corresponding to 
the QuMERA channel $). In other words, exploiting 
the above derivation the constrained energy minimiza- 
tion problem |T]) corresponds to minimize the functional 
F($,p) = Tr[$(p) h] + C \\$(p) - p\\, where we omitted 
the index k as, from now on, all the quantites refers to a 
subsystem composed by three sites and £ is a Lagrange 
multiplier. This is the main result of this paper and as we 
will show in the following, this allows us to compute the 
thermodynamic ground state energy employing a number 
of steps which do not depend upon the system size (the 
limit L — > oo being already included in the QuMERA de- 
scription). Once the transfer operator has been found via 
minimization, as shown in Q, the long range properties 
(critical exponents) can be easily computed. 

A. Minimization 

The problem of finding the ground state properties 
have been reformulated in terms of the minimization of 
the quantity ((3]) with the additional constrain ([2]). To 
solve it we proceed with an iteration procedure starting 
from an initial guess <E>o and Pg for the QuMERA channel 
$ and its associated eigenvector p*. Specifically, all the 
X and A tensors entering in the definition of the trans- 
fer super-operator $ are given by an initial guess xo and 
Ao which we optimize by varying each of them one at 
the time. In this framework the problem is reduced to 
perform the following minimizations 

min{Tr [(AgxAgxAo) 1 Po( a oX^oX^o) h]} , (4) 

where the tensors are contracted following rule given in 
Fig. Q] (here we specialize in the case of the optimiza- 
tion of Xi while keeping fix the initial guess for Ao and 
Pq). Possible strategies to solve such minimization have 
been presented in where unconstrained quadratic op- 
timization have been considered. This approach is not 




FIG. 1: Graphical representation of a one dimensional MERA 
tensor network for a many-body system of L = 8 sites. The 
red elements correspond to the disentaglers tensors \y while 
the blue elements are the isometry tensors A (the green ele- 
ment C is the hat of the MERA which plays no role in the 
L — * oo limit). Any two joined legs from any two distinct 
nodes indicate saturation of the associated indices of the cor- 
responding tensor In the inset it is shown the tensor 
compound which determines the transfer super-operator $ 
associated with the MERA (see Ref. [i| for details). 



longer applicable in our case where the homogeneity con- 
straint adopted to deal with the translational invariancc 
of the system, forces all the xs and As entering in Eq. (g]) 
to be identical while transforming the minimization in a 
non-quadratic optimization problem. To cope with this 
we approach Eq. ^ by computing at each step the lin- 
earized gradient on a basis for the tensors belonging to 
the MERA, while alternating moving from the xs to the 
As. Once the gradient is computed we move the tensor 
in a random direction whose signs in a given basis (going 
forward or backward along a given direction) are defined 
by the gradient's ones as in Ref. (l2l |. 

Once the optimization of the x and the A have been 
performed and a new QuMERA channel <&i is defined, 
the stability constraint ([2]) have to be fullfillcd. Hence, 
we replace pg with a new guess pi defined as the solu- 
tion of the eigenproblem <i>i(p 1 ) = p*, which can be ad- 
dressed by using some smart eigen-problem solvers as the 
Lancsoz or Davidson algorithms. The latter approach is 
favorable also due to the presence of a "good guess" pg 
for the new eigenvector p*. The result presented in the 
next sections are obtained using this second method. We 
can then proceed again in optimizing another tensor be- 
longing to the transfer super-operator up to the desired 
convergence. Notice that this procedure is not garanteed 
to converge to the optimal minimun, however in general 
gives good results. 

III. THE MODELS 

The models we consider are defined through the Hamil- 
tonian 

i 

+ Aa*a* +1 + 2 b a*\ , (5) 
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where af (a — x, y, z) are the Pauli matrices of the «-th 
spin. The constants J, A, 7 and b respectively character- 
ize the interaction strength between neighboring spins, 
the anisotropy parameter along z and in the xy plane, 
and an external transverse magnetic field. The Hamilto- 
nian ((5]) has a very rich structure We consider the 
two cases cases: 

• the Ising-model in a transverse field. Here one has 
A = and 7 = 1. The model presents a critical 
point at 1 6 1 = 1 . 

• the XX Z anisotropic Heisenberg model. Here one 
has b, 7 = and A generic. In this case the Hamil- 
tonian is critical for — 1 < A < 1 while it has fer- 
romagnetic or anti-ferromagnetic order for A > 1 
or A < — 1 respectively. 

The exact critical exponents, related to the correlation 
functions 

K< +e )-K)(< +£ )~r-^ , (6) 

are given by 

• Ising model (7 = 1, A = 0) 

v z = 2, v x = 0.25 , v y = 2.25 , (7) 

• XXZ model (7 = 0, |A| < 1) 

v x = v y = l/v z = 1 — arcos(A)/-7r . (8) 

We are going to compare our numerical results against 
these values. 

IV. RESULTS 

The algorithm we developed for solving the minimiza- 
tion of Eq. ([3]) requires to store order of 0(m 6 ) tensors 
with m being the dimension of the indexes of the x s an d 
As, their moltiplication (0(m 10 ) operations) and their 
diagonalization. As for the t-MERA [lj] this is "reason- 
able" up to m — 4, 6. 

We first concentrate on the Ising model, where we have 
obtained an energy convergence to the exact energy per 
site of the ground state at the thermodynamic limit of 
the order of 5E = 1CP 4 . In Fig. [5] we report (in decreas- 
ing order) the modulus of eigenvalues of the QuMERA 
map $ found by means of the minimization strategy we 
rpesented in the last section. As discussed in Q we ex- 
pect the first eigenvalue (k = 1) to be non-degenerate 
as the $ should be mixing. The subsequent eigenvalues 
instead express the critical exponents of the system via 
the relation: 

\K t a h \=2-"*/ 2 . (9) 

The dashed red lines of the Figure report the theoretical 
expectations of the eigenvalues. The result is good, the 




Progressive Number of the Eigenvali 



FIG. 2: Modulus of the eigenvalues of $ for the Ising model for 
b = 1, J = 1. The dashed lines correspond to the theoretical 
values obtained via Eq. ©. 
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FIG. 3: Energy error (Part A) and critical exponents (Part 
B) of the XXZ model for various values of A. For comparison 
in the inset we report also the other eigenvalues of the map. 
(plots obtained for m = 4). 



first one is exact up to the machine precision, while the 
others have numerical errors of the orderd of 5\K a \ ~ 

icr 4 ,icr 2 ,i(r 2 . 

The results obtained for the Ising model are very 
promising, however it is well known that the Ising model 



4 

















2,5 








2,5 
2 










2 


















V 








1,5 










1,5 
1 

0,5 









V 

0,5 





















10 10" 3 


10" 2 


10'* 


10 3 








A) 






E 


B) 



FIG. 4: Convergence of the cirtical exponents for the XXY 
model for A = 0.05 (left) and A = 0.1 (right). Dashed lines 
gives the exact critical exponents v a . 



has a simple spectrum and properties due to the fact that 
it is equivalent to free fermions [l5|. We then afford the 
study of a more complex model, the XXZ model. In 
Fig. [3] (Part A) we show the error of the energy ground 
state at the thermodynamic limit for various values of A 
and m = 4. As expected, the error is increasing with A, 
however it remains of the order of SE = 10 -3 . In Fig. [3] 
(Part B) we show the correspondent critical exponents. 
As it can be seen, the unital eigenvalue is always present 
(at almost machine precision) while the first two degener- 
ate eigenvalues are related to the v x , v y critical exponents 
(black circle and orange star) with errors increasing with 
A reflecting the precision in the energy of the ground 
state. The third critical exponent is also detected with 
even a better precision than the first two for A > 0.1 . In 
the case A = the eigenvalue correspondent to v z does 
not happen to be the third one as a few more are present 



between the "physical" ones. This behavior should be in- 
vestigated in more details (the prefactors could be zero). 
Finally we report in Fig. [4] the convergence of the critical 
exponents for the XXZ model for different values of A as 
a function of the maximum "move" in the minimization 
procedure e which corresponds to the convergence of the 
ground state. As can clearly seen, the behavior is reg- 
ular, converging to a value that is near to the expected 
one. 



V. CONCLUSIONS 

This paper presents a recursive algorithm for optimiz- 
ing the MERA tensor description of the ground state 
I^g) of a critical many-body Hamiltonian in the ther- 
modynamic limit of infinitely many sites. It is based on 
the results of Ref. Q which established a connection be- 
tween I^g) and the transfer super-operator of the MERA 
network that approximates it. 

While completing the writing of the present 
manuscript we became aware of a work by Pfeifer, 
Evenbly and Vidal (arXiv:0810. 0580^1) employing a 
modified MERA tensor network representation seems 
to lead to a more efficient evaluation of the critical 
exponents. 
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